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Abstract. The strength and stability properties of hierarchical load bearing networks and their 
strengthened variants have been discussed in recent work. Here, we study the avalanche time dis- 
tributions on these load bearing networks. The avalanche time distributions of the V— lattice, a 
unique realization of the networks, show power-law behavior when tested with certain fractions of 
its trunk weights. All other avalanche distributions show Gaussian peaked behavior. Thus the V~ 
lattice is the critical case of the network. We discuss the implications of this result. 
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1. Introduction 

It is well known that many important systems such as the World Wide Web, and the Inter- 
net, power grids, cellular networks can be modeled as complex networks llIlQ]- Branch- 
ing hierarchical networks constitute an important class of networks, and have been used as 
models of granular media jstl, river networks as models of the lung inflation process 
, directed percolation processes |7] and other |8j. Branching hierarchical networks 
have also been studied as models of load bearing networks, and strategies to improve their 
strength @] as well as their tolerance to failure ifioll have been studied in detail. In this 
paper, we study avalanche propagation on the load bearing branching hierarchical net- 
work, and show that the distribution on a specific realization of the lattice, known as the 
V— lattice shows power-law behavior. Thus the V— lattice constitutes the critical case 
of the load-bearing networks. We also discuss the relevance of this result in application 
contexts. 

2. Hierarchical models 

The distributions of avalanche time are studied for a 2 — D hierarchical lattice, henceforth 
to be called the original lattice, two of its variants enhanced by capacity enhancing strate- 
gies i^, and also for a specific realization of the original lattice, called the V— lattice. 

The original lattice is, in fact, the special (z(0, 1) case of the Coppersmith model |3] of 
granular media and Scheidegger's river model ||4|], in which, a site in a layer is connected 
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Figure 1. A network of M=8 layers with 8 sites per layers with connection-probability 
p = 1/2. The beaded line is the trunk of maximal cluster. The weight bearing capacity 
of the trunk is Wt=99. 



randomly to one of its two neighbors in the layer below, yielding a hierarchical river like 
branched structure. We discuss here the version of the lattice set up in Ref. [9]. Each site 
in the lattice has the capacity to support unit weight if it is not connected to any site in the 
layer above, and has capacity w + 1 if it is connected to sites whose capacities add up to 
w, in the layer above. Thus, the capacity w{i^'^) of the i*'' site at any layer Af*'* is given 
by 

where i/^^^ and are the left and right neighbors of the site i*^, in the M — 1th 

layer. The quantity ,1^^^) takes the value 1 if a connection exists between 

and i'^ and if a connection does not exist. FiglT] shows a realization of a original 
network with lattice side M — 8. It is clear that the network consists of many clusters, 
where a cluster consists of a set of sites connected with each other The largest such 
cluster of sites in the network is called the maximal cluster of the network. FiglT]shows a 
representation the network where the connections between sites are shown by solid lines, 
and the numbers in the brackets indicate the capacity of each site. Ci , C2, C3 and C4 are 
the clusters seen in this reahzation of network, and C2 is the maximal cluster. The beaded 
line denotes the trunk of maximum cluster, where the trunk is the set of connected sites 
with highest weight bearing capacity in the maximal cluster. 

2. 1 The V— Lattice Network 

The original lattice has a very special realization which bears the maximum trunk capacity 
compared to all possible realizations. This lattice bears a unique V— shaped cluster that 
includes all the sites in the first layer, and {M — / + 1) sites in the /— th layer, where M 
is the total number of layers. One of the arms of the V constitute the trunk, and all other 
connections run parallel to the arm of the V that is opposite to the trunk. Thus, this cluster 
includes the largest number of sites, and is thus the largest possible cluster the original 
lattice could have. We call this lattice the V~ lattice, and the cluster the V~ cluster. The 
T^-lattice realization with lattice side M = 8 is shown in Fig. |2] The largest cluster here 
is the cluster C2. Every site in the cluster at the layer / has capacity /, except for the 
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Figure 2. The V— lattice network of D=S layers with 8 sites per layers which is the 
critical case of our original network. The beaded line is the trunk of maximal cluster 
with weight bearing capacity Wt=120. 



trunk site which has capacity wt{I) = wt{I — 1) + / in the /-th layer. We will see that 
this network corresponds to the critical case of the original lattice. 

3. Probability distributions of avalanche times 

It is interesting to study the avalanche phenomena on the network. Our study is of rel- 
evance to any situation where threshold phenomena propagate on a network. Examples 
of this range from power propagation on grids, electrical impulses on neural networks, 
ventilation in respiratory networks to directed percolation and granular media. 

We define avalanches on the network in terms of weight transmission. The weight 
transmission in the network takes place along the connections between sites. When a site 
in the first layer of the network receives a weight W, it retains an amount equal to its 
capacity Wc and transmits the rest, i.e. W — Wc, to the site it is connected to, in the layer 
below. Thus, the weight transmission is in the downward direction and the sites involved 
in this process with their connections constitute the path of transmission. Let P be one 
such path and Pjj be the site on P in the D^^ layer. Then, the excess weight at a site Pd 
in the D*^ layer is given by: 

D 
K=l 

If W'^^{Pd) < 0, then the transmission ends at the layer of the path P and is 
considered to be successful. On the other hand, if W^'' (Pd ) > 0, the weight is transferred 
to Pd+1- Finally, if there is still excess weight left at the M*'' layer, it is then transmitted 
to the corresponding site in the first layer and the second cycle of downward transmission 
begins as described above. This process continues tiU either there is no excess weight 
left,i.e. the transmission is successful, or the receiving site is not able to transmit the 
excess to the site in the layer below, and the transmission fails. Such a failure occurs 
when the transmitting site is connected to a site that has already received its share of the 
weight (i.e. saturated it's capacity) in the first cycle of transmission, thus making further 
transmissions impossible. This process of weight transmission is defined as an avalanche. 
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The time taken for an avalanche is defined to be the number of layers traversed by the 
weight in the network. 

In this paper we compare the probability distribution of avalanche times t between the 
original lattice, and its strongest realization, the V— lattice. This probability distribution 
as the name suggests is, in fact, the distribution of the number of layers traversed during 
all cycles of successful avalanche transmission by a test weight placed at a random site 
in the first layer for any lattice. This probability distribution in case of the original lattice 
tested with the trunk weights, has been studied by Janaki et. al |@]. 
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Figure 3. Scaled probability distribution of avalanche time t for original lattice corre- 
sponding to 1000 realizations of networks of side Al — 100, and M = 300. Both, 
the distributions for M = 100 and Al — 300 scale by the total number of layers M in 
network. 



The probability distribution of avalanche times t for the original lattice is shown in 
Fig|3]for an ensemble of 1000 successful weight transmissions for weight equal to trunk 
capacity Wt- In the original lattice no avalanche is seen when t/M < 1 as there are 
no paths with capacity greater than trunk capacity Wt i- e. there are no successful trans- 
missions for t/M < 1. Also the avalanches of weight transmission in original lattice 
can cycle as much as thrice through the network. When these distributions for original 
networks of different sizes are scaled by their respective total number of layers M, they 
collapse on one another as shown in Fig. [3] Similar behavior is seen in Ref. 

The distribution of avalanche times t for original lattice for test weights which are 
fractions of trunk capacity Wt, is different from that for weights equal to trunk capacity 
as shown in Fig. |4] Fig. |4] shows the distribution P{t) tested for weights ranging from to 
O.IWt to 0.9Wt- In the distributions corresponding to 0.1T4^t and 0.2Wt as shown in 
Fig. |4](a), there is only one cycle in the form of a Gaussian peak, which can be expressed 
by the equation 

P{t) = ^exp(-^^-^) (2) 

where a is a constant and cr is the standard deviation. The values of tj and chi-squared 
tested for the accuracy of the fits corresponding to O.lWr and 0.2Wt are shown in the 
caption of Fig. |4] In the distributions corresponding to O.SWt and 0.4Wt the emergence 
of second cycle is seen as shown in Fig. |4](b). This second cycle persists till 0.9Wt, 
besides, it gradually attains the form of a new Gaussian peak as can be seen in Fig. |4](c) 
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Figure 4. Probability distributions of avalanche time t for original lattice corre- 
sponding to 1000 realizations of networks of side M = 100, fit to Gaussian given 
by Eq.(2) when tested for weights equal to (a) O.IWt ((t=4.46, x^=0.0028) and 
0.2Wt (cr=6.08, x^=0.0024), (b) O.SWt (cr=6.92, x^=0.0016) and 0.4H^t (cr=7.34, 
=0.0024), (c) O.GWt (cri=7.75, (J2=17.55, x^=0.0013) and 0.8Wt (cri=4.46, 
(72=14.68, x^=0.0062), and (d) 0.9Wt ((t=19.82, x^=0.0035) and Wt- 



and (d). In the distribution corresponding to 0.9Wt, the first cycle completely loses its 
Gaussian form, and a new third cycle emerges out after the Gaussian peak corresponding 
to the second cycle, which is apparent from Fig. |4] (d). Thus as the test weight start 
approaching the trunk capacity, we start seeing more and more transmissions which cycle 
more than once through the lattice. 

The distribution of avalanche times for the V— lattice is found to be quite different from 
that of the original lattice. It does not show any systematic behavior for weights equal to 
trunk capacity Wt i- e. distinct lattice realizations show distinct behavior However, the 
V— lattice does show interesting behavior for avalanche time distributions when it is 
tested for weights less than its trunk capacity. Fig. |5]displays the probability distribution 
of avalanche times for the V— lattice (1000 realizations). It is clear from the figure that 
the V— lattice displays a power law behavior of the form P{t) ^ with exponent 
a = 2.45 and a = 2.96 when weights equal to O.lWr (Fig. |5](a)) and 0.2Wt (Fig. |5] 
(b)) are placed in the first layer, respectively. The power law regime in the distribution of 
avalanche times gradually starts disappearing when the distribution is tested with values 
of weights placed approaching trunk capacity, and at the trunk capacity no stable behavior 
in the distribution for different realizations is seen. The difference in the behavior of the 
distribution of the V— lattice from the original lattice is due to the presence of a unique 
and asymmetric cluster, called the V~ cluster, in the V~ lattice network. The V— cluster 
includes sites of all levels of capacity. Hence transmissions on the lattice can achieve 
success at any one of the layers. This accounts for the power law distribution and thus 
brings criticality to the distribution of avalanche times t. When the distribution for the V— 
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lattice is tested for weights higher than 0.2Wt the power law regime of the avalanche time 
distribution starts dying down slowly much before the instability in distribution is seen at 
trunk capacity Wt- This behavior of the existence and subsequent disappearance of the 
power law regime in the distribution is one of the indications that the V— lattice is indeed 
a critical case of the original lattice. Such behavior has not been seen in any of realization 
of the original lattice except the V— lattice. Even two connectivity strategies studied by 
Janaki et al. ^ does not show the existence of such behavior. 

The study of avalanche times t of successful weight transmissions along the trunk path 
in the V~ lattice is shown in Fig. |6] From this figure it is clear that the avalanche time 
t along the trunk is also governed by a power law t ^ Wj with /3 = 0.33, against Wf 
which are test weights ranging from 0.1 Wr to Wt- 




Figures. Probability distributions of avalanche time t corresponding to 1000 realiza- 
tions for the V— lattice network of side M = 100 when tested for weights equal to 
(a) O.IWt, and (b) 0.2Wt- Small regimes for O.lWr (as shown in inset of (a)) and 
0.2Wt (as shown in inset of (b)) display power law behavior with exponent a = 2.45 
and = 21.136, and a = 2.96 and = 52.197 respectively. 




log(Wf) 

Figure 6. Avalanche time t displays power law behavior t ~ Wf^ with /3 — 0.33 and 
= 0.00001 against weights placed as percentage increase in trunk capacity for V~ 
lattice for M = 100. 
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4. Conclusions 

Threshold phenomena propagating on branching hierarchical lattices are of interest in a 
diverse variety of application contexts. Models of this kind have been proposed in the 
context of respiration networks jslSlj voter models |8], granular media [3], power net- 
works ifTHl . river networks 0], as well as directed percolation fj^ contexts. Avalanche 
distributions in these contexts have significance for ventilation strategies for respiratory 
networks, opinions or preferences cascading on voter networks as well as numerous per- 
colation contexts. It is interesting to note the significance of the V— lattice configuration 
as the critical configuration in these contexts. It is clear from our results that a V— shaped 
structure which spans the breadth of the lattice at one extreme, with the trunk, i.e. the col- 
lection of the strongest sites along one arm, is best able to support cascades of all possible 
scales. The utility of this configuration for optimizing the connectivity structure of net- 
works deserves to be explored further. We hope to explore some of these directions in 
future work. 
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